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The Solar System bounds on -R" gravity are often ignored in the literature by invoking the 
chameleon mechanism. We show that in order for the latter to work, the exponent n must be 
ridicolously close to unity and, therefore, these theories are severely constrained. 
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The acceleration of the cosmic expansion discovered 
with type la supernovae 1] still lacks a satisfactory expla- 
nation. The hypothetical dark energy wich is supposed to 
drive this acceleration is an ad hoc explanation: it cannot 
be detected directly in the laboratory and is extremely 
exotic due to its negative pressure P. Its equation of state 
should be P ~ — p (where p is the comoving energy den- 
sity) and phantom energy, which opens the door to much 
trouble with its instabilities and thermodynamical be- 
haviour, is not at all excluded by the observations. Much 
theoretical effort has gone into proposing an abundance 
of models for dark energy and to constrain it observation- 
ally (see [1] for a detailed discussion and for references). 
An alternative approach consists of dispensing with dark 
energy and postulating, instead, that Einstein's theory 
of General Relativity (GR) fails at the largest scales and 
that, with the cosmic acceleration, we have detected de- 
partures from the expected GR behaviour. This proposal 
[3,01 has led to a revival of /(i?) or "modified" gravity, 
described by the action 
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where g is the determinant of the spacetime metric gab, 
R is the Ricci scalar, k = SttG, G is Newton's constant, 
and S"*^™) is the matter action. This class of theories, 
which reduces to GR for a linear function f{R), comes 
in three versions: metric, Palatini, and metric-affine for- 
malisms (see [1,0 for reviews and f?\ for introductions). 
The more complicated metric-affine formalism 8^ is not 
fully developed yet and has seen little use in cosmology. 
Inside matter the Palatini formalism, in which the metric 
and the connection are treated as independent variables, 
is riddled with problems unless its field equations get 
modified by higher order terms [9| and, therefore, we will 
discuss here only the metric formalism, in which the con- 
nection is the metric connection (the distinction between 
metric and Palatini formalisms is irrelevant for GR, but 
the two variations produce inequivalent field equations 
for non-linear f{R) functions). 

Metric f{R) gravity contains a scalar degree of free- 
dom, identified with = f'{RX In fact, metric f{R) 
gravity is a Brans-Dicke theory [lOj with parameter a; = 
and a special potential for the Brans-Dicke field </> [ll|. 



Starting from the action ([T]) and introducing a new field 
X, the action 

5=^1 d^xV^ [fix) + fix) {R - X)] + ^^"^ (2) 

is dynamically equivalent to Variation with respect 
to X yields fix) (i? - x) = and X = if /"{R) + 0, 
and the action (jT]) is reproduced. If we define the field 
= /'(x) and set 



the action becomes 



S 



V{(j))] + 
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an w = Brans-Dicke theory [10l|. 

Many choices for the function f{R) have appeared in 
the literature, and there are viable ones which satisfy 
both theoretical viability criteria (such as correct cos- 
mological dynamics, smooth transition between different 
cosmological eras, well-posed initial value problem, sta- 
bility, correct weak-field limit and dynamics of cosmo- 
logical perturbations) and experimental constraints 5, 6j . 
There is a large body of literature ([l3,[l3| and references 
therein) on the choice f{R) = aR^ (where a > has the 
dimensions of a mass squared and n is not restricted to 
be an integer), on which we focus. Let us be clear on 
the terminology here: often, the literature refers to the 
theory described by f{R) = R + aR"^ motivated by quan- 
tum corrections to the Einstein-Hilbert Lagrangian as 
"i?^-gravity" (and, consequently, to f{R) = R + aR^ as 
"i?"-gravity" ) . This is not what we mean here: the term 
"i?"-gravity" in this paper refers strictly to the choice 
/(i?) = ai?" and our considerations apply only to this 
class of theories (the prospects appear much better for 
f{R) = R + aR"^ theories). 

i?" gravity, like any f{R) theory, is subject to experi- 
mental constraints: while, from the mathematical physics 
point of view, it is perfectly acceptable to study this the- 
ory as a toy model in order to obtain analytical or qual- 
itative insight on exact solutions, or on the role that the 
scalar degree of freedom f'{R) may play in modifying 
GR, or to replace the full theory f{R) — R + aR" with 
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ai?" (which is mathematicahy easier to handle) for n > 
in the strong gravity regime, the exponent n is not an en- 
tirely free parameter if the theory is meant to constitute 
a realistic alternative to dark energy. 

Let us first consider two basic theoretical requirements 
associated with stability. First, avoiding the Dolgov- 
Kawasaki instability ^ leads to f"{R) > [H,]!?! 
which corresponds to n < or n > 1 for f{R) = aR^ . 
Second, on time scales shorter than the Hubble time, 
one models the present universe as a de Sitter one, and 
de Sitter space is usually found to be a late-time attrac- 
tor in f{R) and dark energy models. Therefore, it is 
important that de Sitter space be stable, too, which pro- 
vides the second criterion. The mass m of the scalar field 
(f) = f'{R) in de Sitter space is given by 
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where a zero subscript denotes quantities evaluated in 
the de Sitter space with Ricci scalar i?o. Eq. ([5]) has 
been derived in a variety of ways, including the weak- 
field limit 
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- 20| , gauge- invariant perturbation analyses 
of de Sitter space [21| , and calculations of the propagator 
of f{R) gravity in a locally flat background [2^. For 

f{R) = aR"^, it is 
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O/ (6) 

and the requirement that the field (j) be non-tachyonic is 
equivalent to 1 < n < 2. We take the parameter n in the 
intersection of these two intervals 1 < n < 2 (bounded 
from below by GR). 

In order for /(i?) — ai?" to provide a realistic alterna- 
tive to dark energy, it also needs to satisfy the available 
experimental constraints. Wri ting n = 1 + S, light deflec- 
tion does not provide bounds [2 



24{ but the precession 



of Mercury's perihelion yields the stringent limits j23l.l25 



(2.7 ±4.5) • 10' 
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This constraint is often ignored in studies of i?" gravity 
[l^ . based on the belief that the Solar System limits are 
circumvented because in the weak-field limit of general 
f{R) gravity, the effective degree of freedom cj) = f'{R) 
is endowed with a range which may be very small at Solar 
system densities and much larger at cosmological densi- 
ties. This feature would enable effects on cosmological 
scales but would shelter (f) from the experimental bounds 
in the Solar System (the chameleon mechanism at work, 
see below). This argument is misleading: let us examine 
how it applies to the weak field limit of f{R) gravity in 
general, and then discuss the specific i?" theory. 

The weak-field limit of f (R) gravity has been studied 
by various authors 0, [s^, [s^l ■ Based on the equiva- 
lence between metric f{R) and uj = Brans-Dicke grav- 
ity and on the Cassini bound > 40000 ^28i], early work 



dismissed all f{R) theories as unviable [29[. However, the 
fact was missed that the Cassini limit only applies to a 
Brans-Dicke field with range larger than, or compara- 
ble to, the size of Solar System experiments, while the 
effective mass and range of the scalar field cj) — f'{R) 
depend on the background curvature R, hence on the en- 
ergy density of the environment. This is the chameleon 
mechanism originally discovered in quintessence models 
of dark energy [s^, and later rediscovered in modified 
gravity [s^ j. The chameleon mechanism is not imposed 
to fine-tune the theory and evade the experimental lim- 
its: it is contained naturally in f{R) gravity and whether 
it works or not depends on the specific theory considered. 

In the weak-field limit of f{R) theories [2^, [soj, one 
considers a spherically symmetric, weakly gravitating, 
perturbation of mass M of a cosmological space. In an 
adiabatic approximation, the background is taken to be a 
de Sitter space (with constant curvature, Rab — RoQab/'^, 
and i?o = 12_ffQ), which is a solution of f{R) gravity sub- 
ject to the conditions j5|] 
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JqRq = 2/o . Hq — A — — 



The weak-field line element is written as 



[1 + 2$(r) + H^r"^] dr^ 



r'^dn^ . 



(9) 



where dfi^ = d9^ + sin 



^ is the line element on the 



unit 2-sphere and ^>{r) and <I>(r) are post-Newtonian po- 
tentials. The goal is to compute these potentials by solv- 
ing the linearized fourth order field equations and to ob- 
tain the PPN parameter 7 = — ^1^/$, which is subject to 
the Cassini bound [281 



I7- 1| < 2.3 • 10" 



(10) 



(in GR, vjf = -$ is the Newtonian potential -KM/(87rr) 
and 7 = 1). A linearized analysis assuming that 
|*(r)| , |$(r)| < 1, Har < 1 f{R) is analytical at i?o, 
and mr < 1 yields [ifl, S, 113 



kM 



$(r) = 



kM 



1=\, (11) 



in gross violation of the Cassini bound. This result would 
spell the end for /(i?) gravity if it wasn't for the fact 
that the assumption of a light scalar field, mr ^ 1, is 
violated. In many f{R) theories this happens naturally 
and the mass of is large at high (i.e.. Solar System) 
densities and almost zero at cosmological densities [34| . 
But does this mechanism work for /(i?) = aR"7 

To answer this question, note that for n — 2 (the 
largest value of n allowed by theoretical stability) the 
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mass dni) of the scalar (j) = f'{R) in a de Sitter back- 
ground is exactly zero and this field has infinite range in- 
dependent of the density of the environment; therefore, 
it is certainly subject to Solar System constraints and B? 
gravity is ruled out experimentally. 

At a first sight, it looks surprising that the mass m van- 
ishes while the potential ^ turns out to be V{4)) — 
for this theory. The solution to this apparent contradic- 
tion is that it is not ((/>), but rather the combination 
(j)^ — 2V{<j)) that enters the equation of motion for the 
Brans-Dicke scalar [10| 
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where T^™) is the trace of the matter stress-energy tensor 
-^i™^ (which, in the weak-field, slow-motion limit, reduces 
to —p) and (f>^ — 2^(0) vanishes identically for a purely 
quadratic potential (isj . 

Incidentally, the theory f{R) = aR^ with a > (in 
D spacetime dimensions, /(i?) = ai?^/^ J[36|]) has other 
peculiarities or theoretical problen is ISTI . ISSj : it does not 
have the correct Newtonian limit 39] and eq. ^ is sat- 
isfied for all, not for special, values of the Ricci cur- 
vature R, which leads to unpleasant consequences (40| . 
That something goes wrong in the weak-field limit can 
be seen in the post-Newtonian potentials dTTI) which, us- 
ing Ro 



127?Q, reduce to 



^ = -2$ = - 
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1 R,cH^' cHn 
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(13) 

(restoring G and c) where Rg = 2GM / is the 
Schwarzschild radius of the mass M. 4' and $ are no 
longer guaranteed to be small in absolute value because 
cH^^ /r <^ 1 and it is not clear how to choose the param- 
eter a. A more refined analysis including terms of order 
Hnf yields post- Newtonian potentials with Yukawa terms 



* = -- 



GM 



1 - 



5 c- 



GM ( 5(l + ar)c- 



(14) 
(15) 



In the limit n — > 2^ in which a — ?> and the range of 
the scalar becomes infinite, the Yukawa terms dominate 
the Newtonian ones and diverge. The range of must be 
kep t small in order to recover even the Newtonian limit 
44|. 

At the opposite range of values for n we have GR, 
which is viable and in agreement with all available Solar 
System experiments. Between the values n — \ and n — 
2, the range of the scalar field varies continuously but 
rapidly between zero and infinity. This range is given by 



0.75 




0.25 



FIG. 1: The range s(n) of the effective scalar degree of free- 
dom ^ = f{R) (in units cH^^) for the theory f(R) = ai?". 
The parameter n is in the range 1 < n < 2 allowed by sta- 
bility. The function s(n) starts out with vertical tangent at 
n = 1. 



the function 



sin] = — ,, 
^ ' 2 V 2-n 



(16) 



in the interval [1,2]. This function varies continuously be- 
tween s(l) = and its limit lim5_^2- ■5("-) = +oo, always 

increasing. The derivative s' [n) = 2{n-\)^\'2-nfl'^ 
ways positive and the tangent to the graph of s(ri) starts 
vertically at n = 1 and ends vertically as n — > 2", which 
means that the range of increases quickly as the i?"- 
theory departs very slightly from GR (see fig. [Ij . Clearly, 
as long as the exponent n is very close to unity, the theory 
behaves as GR and passes the experimental tests while, 
approaching values of n closer to 2, the experimental 
bounds begin being violated, and disaster happens in the 
limit n — )• 2^. In conclusion, only for very small values of 
n it is possible to invoke the chameleon mechanism in the 
weak-field analysis. By imposing the range of the scalar 
to be less than one astronomical unit (1.496 • 10^^ cm) 
and using the value TJq = 70 km • s~^ • Mpc~^ for the 
Hubble parameter, one would obtain the requirement 



< (5; 



1< 5 • 10 
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(17) 



Of course, realistic Solar System experiments do not have 
this level of precision, and the limit ([7]) applies instead. 
This renders i?" gravity a poor candidate for a realistic 
alternative to dark energy. 
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